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Abstract 

We study the stabilization of coherent suppression of tunnehng in a driven double-well system 
subject to random periodic (^—function "kicks". We model dissipation due to this stochastic process 
as a phase diffusion process for an effective two-level system and derive a corresponding set of Bloch 
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' equations with phase damping terms that agree with the periodically kicked system at discrete 

o ■ 

. times. We demonstrate that the ability of noise to localize the system on either side of the double- 



o 



well potenital arises from overdamping of the phase of oscillation and not from any cooperative 
effect between the noise and the driving field. The model is investigated with a square wave 



■ drive, which has qualitatively similar features to the widely studied cosinusoidal drive, but has the 

^3 : 

Q^, additional advantage of allowing one to derive exact analytic expressions. 
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I. INTRODUCTION 



Coherent suppression of tunneling is a localization effect that occurs when a potential with 

multiple minima is exposed to a periodic drive in a specific parameter regime [0, 0, 0] . In their 

original paper, Grossman et. al. 0] examined this phenomena in the context of a quartic 

double-well potential driven by a strong cosinusoidal force. It was later shown that the 

essence of this effect could be described by the two lowest energy eigenstates corresponding 

to the symmetric IS"), and antisymmetric 1^4), states of the ground state doublet. This 

approximation is valid if the drive strength, drive frequency and energy splitting are small 

compared to the energy spacing between the average energy of the ground doublet and the 

higher energy levels Q . The two level approximate Hamiltonian for this system written in a 

basis of left and right states, \L,R) = {\S)^\A)) / \/2, is described by a bare tunneling system 

with an energy splitting and a cosine driving force with amplitude and frequency uod, 

A A 

H{t) = -Y^x + cos{uJdt)a^. (1) 

A complete suppression of tunneling occurs when the drive amplitude is much larger than 
the energy splitting {A^ ^ A^) and the ratio of the drive amplitude to the drive frequency 
is equal to a root of the zero order Bessel function 0]. This absence of coherent oscillation 
is due to the two Floquet states of the driven system becoming exactly degenerate 0. 

A counterintuitive effect P, H, 0] occurs when dissipation is added to this system. This 
can be modeled by the addition of a stochastic 5-function kick, which is applied periodically 
0. When the ratio of the amplitude to the frequency of the periodic driving field is near but 
not equal to the Bessel root condition in the absence of noise, coherent oscillations between 
the left and right wells proceed at a modified tunneling rate corresponding to the quasienergy 
splitting (See Fig. |T]). At first, the effect of noise is to destroy coherent oscillations causing 
the system to evolve into a 50/50 mixture of left and right wells. However, increasing 
the noise strength appears to have a stabilizing effect on coherent suppression of tunneling 
IP, P] causing the system to remain in the well corresponding to its initial state. This 
seems surprising because the condition for coherent suppression of tunneling depends on the 
state of the system being in a stationary superposition of the two degenerate Floquet states 
for all times whereas dissipation generally tends to destroy quantum coherence. However, 
we will demostrate that this localization effect can be understood intuitively as arising 
from strong damping of the phase of oscillation. This has been previously pointed out by 
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Makarov [§] using a semiclassical quantization scheme to establish a correspondence between 
a classically driven two-level system and the spin-boson problem. Standard treatments of 
dissipation based on the master equation of a reduced density operator after tracing over a 
large reservoir ^ |^ can obscure the physical content of this localization phenomena. Our 
treatment uses a simplified model based on a damping due to phase diffusion from random 
classical kicks. This model corresponds to the regime of weak coupling to a classical noise 
reservoir P, |lOl. From this we derive a detailed and exact expression for the tunneling 



oscillations that gives a clear interpretation of this phenomena that excludes any cooperative 
effects between the noise and the driving field. 

Dissipation in driven two level systems has been extensively studied 0, ^, H, ^ [l^, 



13| , p^ , p!5| , |16|, |T^, |1^, |T9[ because this model provides a realistic paradigm for understanding 
a variety of physical systems. For example this model can be used to describe charge oscilla- 



tions in semiconductor double wells pO, BJ, 22 1, magnetic flux dynamics in superconducting 



quantum interference devices (SQUIDs) [Q, condensed phase electron and proton reactions 



25| , |26| and strong field spectroscopy Recently, quantum tunneling over meso- 



scopic distances and the generation of mesoscopic superposition states has been realized in 
an optical lattice of double wells . This system is ideally suited for the study of the driven 
double well in the presence of dissipation given the large degree of experimental control one 
has over the system For example the energy barrier and energy assymmetry can be 

dynamically controlled through laser beam configuration and externally applied magnetic 
fields. Furthermore, this system can operate in an essentially dissipation free environment 
when the lattice lasers are sufficiently far detuned from the atomic resonance. Dissipation 
can then be reintroduced into the system in the form of well-controlled fluctuations in the 



potential 0, ^ 



In this paper we replace the cosine drive field with a periodic square wave drive field. 
Our Hamiltonian reads, 



Hit) = + ^A{ujdt)a, 



(2) 



J+1 for < mod.,(t) < f 
[ -1 for f < mod^,(t) < r^, 

where = ^ is the drive period. This drive has similar qualitative features to the sinusoidal 
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force. For instance, the condition for coherent suppression of tunnehng occurs when the ratio 
of the amphtude to the frequency of the drive is equal to an integer rather than a root of the 
zero order Bessel function. However, this Hamihonian has the advantage that one can derive 
analytic results without making approximations that are common when analyzing sinusoidal 
drives. The paper is organized as follows. In Sec. II we present an effective magnetic field 
formalism in the Floquet basis showing how one can describe the tunneling system with a 
set of Bloch equations associated with spin precession about a fictitous static magnetic field. 
In Sec. Ill we show how one can use a phase diffusion model to describe the dissipation in 
terms of dephasing terms added to the Bloch equations and we present results that compare 
numerical simulations to analytic expressions. In Sec. IV we summarize our results. 



II. THE EFFECTIVE MAGNETIC FIELD IN THE FLOQUET BASIS 

As discussed above, we simplify the problem of a tunneling wave packet in a double well 
by considering the problem restricted to the two lowest energy eigenstates. Because this is 
a two-level system one can visualize the dynamics geometrically on the Bloch sphere with 
spin up and spin down along the quantization axis z corresponding to the localized \L) 
and l-R) states respectively. According to the bare Hamiltonian, the first term in Eq. (2), 
tunneling on the Bloch sphere is pictured as the Larmour precession of a Bloch vector, s, 
about a static magnetic field in the x direction with a frequency ujl = {h = 1). When 
a drive is added to the tunneling system, such as the square wave drive given in Eq. (2), 
the Hamiltonian becomes time dependent. In general, a drive term causes trajectories to 
explore the entire Bloch sphere in a complicated fashion (See Fig. |^). However, when the 
ratio of the amplitude to the frequency of the drive is an integer and ^ A,j., one finds 
that the trajectories on the Bloch sphere form closed "loops" near the top (or bottom) of 
the Bloch sphere which correspond to the quantum system remaining localized in the \L) 
(or \R)) state (See Fig. |^). This is the coherent suppression of tunneling. Note that if the 
condition A^ ^ A^ (typically A^/A^ ^ 5) is not satisfied then the system can still form 
closed trajectories but they do not stay "near" the top or bottom of the Bloch sphere and 
hence are not localizing. 



According to the Floquet theorem |32, 33], a system with a Hamiltonian that is peri- 



odic in time, H{t + r^) = H{t), has solutions to the Schrodinger equation, \ej), which are 
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eigenfunctions of the single period time propagator or Floquet operator, = 1) 




U{Td,0)\ej) = exp{-iejTd)\ej), 



(4) 



where ej is a quasienergy which belongs to a family of quasienergies such that Sj + kud (where 
k is an integer) belongs to the same physical state. Quasienergies can be uniquely defined 
by requiring them to continuously approach the energies of the time independent Hamil- 
tonian as the periodic part vanishes. Because the Floquet states are stationary states at 
integer multiples of r^, we restrict time to these discrete values. In this stroboscopic picture 
of the dynamics, the Floquet states evolve like energy eigenstates for a time-independent 
Hamiltonian with the quasienergies playing the role of energies. 

For a two-level system we denote the Floquet states Writing the single 

period time propagator in the Floquet basis and neglecting global phase factors gives. 



where Ae = £+ — e_ is the quasienergy splitting and the quantization direction, ep, is chosen 
along \e+). We see that the application of the single period time propagator can be thought 
of as a rotation about an effective magnetic field, = Ae ep. Because the single period 
time propagator commutes with itself at integer multiples of the drive period, the unitary 
time propagator after £ drive periods is found from £ applications of the single period time 
propagator, U{£Td, 0) = f/(rrf, 0)^ (£ = 0, 1, 2, ■ ■ ■ ). The end result of utilizing this represen- 
tation is that one can replace the continuous dynamics of a time dependent Hamiltonian 
with a discrete dynamics associated with a time independent Hamiltonian described by spin 
precession about a fictitous static magnetic field (See Fig|^). This transformation is not an 
approximation, like the well known rotating wave approximation; the two pictures are in 
exact agreement at these discrete times. However, the time resolution of the Bloch trajec- 
tories is fundamentally limited by the drive period because this quasistatic picture does not 
contain information about time scales shorter than the drive period. 

For the specific case of the square wave drive the time dependent Hamiltonian can be 
decomposed into two time independent Hamiltonians, H± = ± ^^z- According to Eq. 
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(2), the single period unitary operator can be easily contructed by multiplying the two half 
period unitaries, 

U{t,,0) = exp(-zi^_^)exp(-zi^+^) (6a) 
= (^l-2sin'(^)sin2(e))i 
-(2i sin(^) cos(^) sin(6')j(T^ 

+ (^2isin2(^)sin(^)cos(e))aj^, (6b) 

where = + and tan(^) = A^/Ax- The effective magnetic field for the square 

wave drive can be found analytically by comparing Eq. (6b) to 

exp(^-i(a-fi)^|^) = icos(^^|^) - z(ct • fi) sin(^^^), (7) 

and solving for the four quantities {Ae, n^., rij^, ra^}. With a proper coordinate transforma- 
tion, we find that the effective magnetic field for a square wave drive can always be viewed 
to be in the x direction, {rix = l,ny = 0,n^ = 0}, and the quasienergy splitting is given by 

Ae = — cos-^ (l-2 sin2(0/2) sm\9)] . (8) 

One can verify that the effective magnetic field vanishes (Bloch vector trajectories form 
closed "loops" on the Bloch sphere) when the second term in Eq. (8) vanishes (sin(0/2) = O) 
giving the condition a/^^ + Al/ujd = j (j is an integer). When A^ ^ A^ this reduces to 
the condition, Az/uJd = j, the same condition stated in Sec. I for coherent suppression of 
tunneling by a square wave drive. 



III. DISSIPATION 

We next add dissipation to the system with stochastic periodic 5- function kicks . Equa- 
tion (2) is modified to read. 

Hit) = + (^^Koodt) + - ^^d))^., (9) 

where is a discrete random variable governed by a distribution, P{C)i defined by a 
mean, (^), and variance, (^^) = cr^. This model can be related to the widely studied 
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spin-boson model in the weak coupling limit where dissipation is accomplished through cou- 
pling to a reservoir of harmonic oscillators at temperature T. We identify the variance with 
(T^ = Ah^^J{Ae)coth.{hAe/2kbT), where J{u!) is the "spectral function" of the reservoir 



Without loss of generality we take {^) = since a non-zero mean can be removed by an 
appropriate redefinition of the single period time propagator 

f/(r,,0)^exp(-2(^-z)^) f/(r,,0). (10) 

The average dynamics are essentially independent of the specific shape of the distribution, 
P(^), because of the central limit theorem. On the Bloch sphere, the dynamics can be viewed 
as single period time rotation about the effective magnetic field followed by a stochastic 
rotation about the z axis by an angle C,e, 

m = (n,{0^np{Aer,)y s(0). (11) 

Taking the rotations about z is completely general; kicks about other directions can be 
implemented with an approprate unitary transformation. The mean trajectory is computed 
by averaging over many realizations of the stochastic evolution. To see how the averaging 
effects the dynamics, consider the average x and y components of the Bloch vector, {sj(t)) 
{j = x,y). After i periods of the drive, the random rotations about z cause the components 
to execute a random walk in the x — y plane. Accordingly, the angle in the x — y plane obeys 
a diffusion equation and the mean x and y components decay as {sj{£Td)) = exp(— ^), where 
the characteristic dephasing time is T2 = ||3^. These decay terms due to dephasing 

can be included in the Bloch equations, 

Sx (^^y ^yj 

Sy = {-Bx s.^ + B:, - (12) 

(^Bx Sy By SxJ ■ 

Even though these are differential equations associated with a continuous time variable, 
there is an implicit understanding that the solution to these equations agree with the real 
system only at multiple integers of the drive period, t = ir^. 

As an example, we show how these Bloch equations can be used to derive the localization 
effect of dissipation caused by stochastic rotations about z that was discussed in Sec. I. We 



take the initial state to be localized in the left well. Because the effective magnetic field for 
the periodic square drive is only in the x direction, the dynamics take place completely in 
the y — z plane. This reduces the number of differential equations that need to be solved to 
two which can be written as a single second order differntial equation, 

+ ^s, + (Ae)^ s,^0 (13) 

with the initial conditions, {^^(O) = 1, Sz{0) — 0}. The solution to this equation is the 
familiar damped harmonic oscillator which can be solved in three regions: underdamped 
(As^ > critically damped (As^ = and overdamped (As^ < where P = l/2r2, 

' ^e-/^* cos(VA£2-/?2 t-5), Ae'> (3' 

(14) 

where tan(5) = f3/^Ae^-(3^. 

For a particle intially in the \L) state and for parameters of the Hamiltonian, — 
0.1, Az — 1.1, and a;^ = 1, we perform numerical simulations according to the procedure 
described by Eq. (11) and average over 500 Bloch trajectories. We compare plots of the 
probability to be in the \L) state vs. time for these simulations with results obtained 
from Eq. (14) and find very good agreement (Fig. 4 a-c). However, the model eventually 
breaks down when T2 < (Fig. 4 d). In this regime the system is decaying faster than 
the oscillation frequency of the drive, the fundamental time resolution for working in the 
Floquet basis as discussed in Sec. II. 

We are now in a position to explain the localization effect of dissipation that was presented 
in Sec. I. For small a, dissipation causes the quantum oscillations to decay on a time scale 72- 
Increasing a (decreasing T2) causes the oscillations to decay on a shorter time scale until they 
reach a critical damping where the coherent oscillations are completely destroyed. Increasing 
a further causes the system to be overdamped with a dephasing time given approximately 
by 1/(T2A£^). In this regime increasing a (decreasing T2) causes the system to decay on a 
longer time scale. The classical noise locahzes the state by destroying the phase coherence 
necessary for transitions between the left and right wells and not by any cooperative effect 
between the noise and the driving field. This can be seen from the fact that this localization 



effect would still occur, for a small enough bare tunneling splitting, even if the drive were 
completely eliminated! The reason this localization effect is stronger "near" the condition 
for coherent suppression of tunneling is because the dephasing time for this overdamping is 
longer at smaller quasienergy splittings. This explains why dissipation appears to "stabilize" 
coherent suppression of tunneling. We have performed simulations for the two-level system 
driven by a cosine drive and have found similar results, the only difference being that one 
must compute the effective magnetic field numerically. 



IV. CONCLUSION 



In this article we studied the localization effect of dissipation caused by random periodic 
5-function "kicks" on a driven two-level system. Utilizing a Floquet formalism we found 
that this type of dissipation could be modeled with a set of Bloch equations including phase 
damping terms. For dephasing times longer than the drive period, these equations yield 
dynamics that agree with the periodically kicked system at discrete times. We considered a 
square wave drive in this article because it gives the same qualtative results as the widely 
studied cosine drive, but it is simpler to study. We found that the tunneling oscillations obey 
the equations of a damped harmonic oscillator and that this localization effect corresponds 
to overdamped tunneling oscillations. Thus, this stabilization of coherent suppression of 
tunneling is due simply to strong phase damping of tunneling oscillations, which effectively 
projects the system into the pointer-basis set by the interaction of the system with the noisy 
environment ||35[ , rather than a cooperative effect between the noise and the driven system. 
The simplified model studied in this article allowed one to gain physical intuition about the 
stabilization of coherent suppression of tunneling in the presense of noise and may prove 
useful in the study of other driven dissipative quantum systems in the weak coupling limit 
such as dynamical localization in a lattice |^, ^ or quantum stochastic resonance ^ 



The utility of the square wave drive over the commonly studied sinusoidal drives may find 
wider application in more general treatments of driven quantum dissipative systems. 
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Fig. 1: The probability to be in the left well, Pl, vs. time, t/r^. For the central cosinu- 
soidal curve (solid line) cr = 0, where a is the root mean variance of the probability distribu- 
tion which governs the discrete random variable, we have coherent oscillations between the 
\L) and \R) states at a modified tunneling rate corresponding to the quasienergy splitting. 
In the central rapidly decaying curve (dashed), a = 0.25, stochastic (5-function kicks rapidly 
destroy coherent oscillations. In the upper, slowly decaying curve (dash-dotted), cr = 2.5, 
the coherent suppression of tunneling has been partially restored by the increased noise 
strength. Compare this to Fig. 7 b on pg. 243 in Ref. ^ and Fig. 2 on pg. 9 in Ref. 

Fig. 2: The representation of coherent suppression of tunneling on the Bloch sphere for 
a two- level system initally in the \L) state: a) Trajectories off the condition of coherent 
suppression of tunneling wind up and down the Bloch sphere, b) When the condition for 
coherent suppression of tunneling is met trajectories form closed "loops" near the top {\L)) 
of the Bloch sphere. 

Fig. 3: Analyzing the system in a Floquet basis allows one to replace the (a) complicated 
continuous dynamics associated with a time dependent Hamiltonian with a (b) simple dis- 
crete dynamics associated with Larmour precession about a fictitous static magnetic field, 
B^. The two pictures agree exactly at these discrete times corresponding to integer multiples 
of the drive period. 

Fig. 4: Comparison of numerical simulations (dotted line) to analytic expression (solid 
line). Each graph plots the probability Pl to be in the left well vs. time, t/r^, for the 
parameters = 0.1, Az = 1.1, and Ud = 1. There is excellent agreement for cases: a) 
underdamped, (cr = 0.20, T2 = ^Ord), b) critically damped (cr = 0.68, T2 = 4.28r(i), and c) 
overdamped (cr = 1.41, T2 = Ir^^). The model breaks down: d) (cr = 2.5, T2 = .32r^) when 

T2 < Td. 
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FIG. 2: Fig. 2, J. Grondalski, Phys. Rev. E 
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